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Abstract 

Moving highly-charged ions carry strong electromagnetic fields that act as a field of photons. In 
collisions at large impact parameters, hadronic interactions are not possible, and the ions inter- 
act through photon-ion and photon-photon collisions known as ultra-peripheral collisions (UPC). 
Hadron colliders like the Relativistic Heavy Ion Collider (RHIC), the Tevatron and the Large 
Hadron Collider (LHC) produce photonuclear and two-photon interactions at luminosities and 
energies beyond that accessible elsewhere; the LHC will reach a jp energy ten times that of 
the Hadron-Electron Ring Accelerator (HERA). Reactions as diverse as the production of anti- 
hydrogen, photoproduction of the p°, transmutation of lead into bismuth and excitation of collec- 
tive nuclear resonances have already been studied. At the LHC, UPCs can study many types of 
'new physics.' 
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I. INTRODUCTION 



In 1924, Enrico Fermi, then 23 years old, submitted a paper "On the Theory of Collisions 
Between Atoms and Elastically Charged Particles" to Zeitschrift fiir Physik lj. This paper 
does not appear in his "Collected Works". Nevertheless, it is said that this was one of 
Fermi's favorite ideas and that he often used it later in life j^J]. In this publication, Fermi 
devised a method known as the equivalent (or virtual) photon method, where he treated the 
electromagnetic fields of a charged particle as a flux of virtual photons. Ten years later, 
Weiszsacker and Williams extended this approach to include ultra-relativistic particles, and 
the method is often known as the Weizsacker- Williams method . 

A fast-moving charged particle has electric field vectors pointing radially outward and 
magnetic fields circling it. The field at a point some distance away from the trajectory of 
the particle resembles that of a real photon. Thus, Fermi replaced the electromagnetic fields 
from a fast particle with an equivalent flux of photons. The number of photons with energy 
u, n(u), is given by the Fourier transform of the time-dependent electromagnetic field. The 
virtual photon approach used in quantum electrodynamics (QED) to describe, e.g. atomic 
ionization or nuclear excitation by a charged particle can be simply described using Fermi's 
approach. 

When two nuclei collide, two types of electromagnetic processes can occur. A photon 
from one ion can strike the other, Figure ^a), or, photons from each nucleus can collide, in 
a photon-photon collision, as in Figure ^b). 

Ultra-peripheral hadron-hadron collisions will provide unique opportunities for studying 
electromagnetic processes. At the LHC, photon-proton collisions will occur at center of mass 
energies an order of magnitude higher than are available at existing accelerators, and photon- 
heavy ion collisions will reach 30 times the energies available at fixed target accelerators. 
The electromagnetic fields of heavy-ions are very strong, so reactions involving multi-photon 
excitations can be studied. 

Ultra-relativistic heavy-ion interactions have been used to study nuclear photoexcitation 
(e.g. to a Giant Dipole Resonance), and photoproduction of hadrons. Coulomb excitation 
is a traditional tool in low energy nuclear physics. The strong electromagnetic fields from 
a heavy ion allow for the study of multi-photon excitation of nuclear targets. This allows 
the study of high-lying states in nuclei, e.g. the double-giant resonance (J, E|. Multiple, 
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FIG. 1: (a) One-photon and (b) two-photon processes in heavy ion collisions, (c) Geometrical 
representation of the photon fluxes at a point outside nuclei 1 and 2, in a collision with impact 
parameter b. The electric field of the photons at that point are also shown, (d) Feynman diagram 
for qq production through photon-gluon fusion to leading order. (e,f) Example of higher order 
corrections to pair-production: (e) Coulomb distortion, and (f) production of multiple pairs, (g) 
The dominant diagram for Au + Au — » Au* + Au* + p° and (h) for Au + Au — ► Au* + Au* + e + e~ 
or a meson X. The dotted lines in panels (g) and (h) show how the mutual Coulomb nuclear 
excitation factorizes from the particle production. 



independent interactions among a single ion pair are also possible. Reactions like multiple 
vector meson production can be used for studies involving polarized photons. The high 
photon energies can be used to study the gluon density in heavy nuclei p at low Feynman- 
x. 

The cross section for photoproduction is 



ox 



n (u>) 



.1 
X 



(1) 



where a\ (u) is the photonuclear cross section. 
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Photon-photon (or "two-photon") processes have long been studied at e + e~ colliders. 
They are an excellent tool for many aspects of meson spectroscopy and tests of QED. 
At hadron colliders, they are also used to study atomic physics processes, often involving 
electrodynamics i n stro ng fiel ds. O ne striking success was the production of antihydrogen 



atoms at CERN's 



17l| LEARI72I Q and at the Fermilab Tevatron [8J. At the highest 



energy colliders, reactions like 77 — > X may be used to probe the quark content and spin 
structure of mesons resonances. Production of meson or baryon pairs can also probe the 
internal structure of hadrons. At the LHC, electroweak processes such as 77 — > W + W~ 
may be probed. The cross section for two-photon processes is ^ 

[a a n{u x )n{u 2 ) 77 , , , , 

J &i Ll>2 

where a^ 7 (ui,u 2 ) is the two-photon cross section. 

Fermi's method has found application beyond the realms of QED. It has been extended 
to strong interactions between nuclei in peripheral collisions. These interactions are mainly 
mediated by pion exchange, and an equivalent pion method has been applied, to describe 
subthresll0ld pion production in nucleus-nucleus colons Q. Fe^acK used the te™ 
nuclear Weiszsacker-Williams method to describe excitation processes induced by the nu- 



clear interaction in peripheral collisions of heavy ions 



111 ]. More recently, a non-Abelian 



Weiszsdcker- Williams field was used to describe the boosted gluon distribution functions in 



nuclear collisions 
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Since Fermi's original work, much progress has been achieved in this field, especially with 
the advent of relativistic heavy ion accelerators like the Bevalac accelerator at Lawrence 
Berkeley National Laboratory (LBNL). Intermediate energy proc esses have been explored at 
heavy ion accelerators at NSCL/MSU, GANIL, RIKEN, and GSl(l7aj. These facilities have 
explored the collective excitation and electromagnetic fragmentation of nuclei, and studied 
many reactions that occur in the sun, supernovae, and the big bang. Experimental studies 
of higher-energy processes have recently begun at Brookhaven's RHIC. These studies have 
included vector meson spectroscopy and production of e + e~ pairs. In the next few years, 
CERNs LHC will begin operations, allowing for the study of heavy mesons, measurements 
of gluon distributions in nuclei, and searches for a host of 'new physics' processes. 

This review will discuss these experiments, their theoretical interpretation, and some 
future possibilities in this field. UPCs have been previously reviewed by a number of authors 
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FIG. 2: Highly energetic charged particles have Lorentz contracted electric fields. The interaction 
of these fields can be replaced by the interaction of real (or quasi-real) photons. 




A. The Photon Flux 



The flux of equivalent photons from a charged particle is determined from the Fourier 
transform of the electromagnetic field of the moving charge. The fields of a relativistic 
particle Lorentz contract toward a co-moving pancake, as is shown in Fig. El The photon 
energy spectrum depends on the time a target particle spends in this pancake, i.e. on 
the minimum distance between the target and the charge and on the projectile velocity; 
the minimum photon wavelength is the width of the pancake at the target. At an ion- ion 
separation (impact parameter) b, the interaction time is At ~ 6/(71;). In the lab frame, the 
maximum photon energy is 

U = At ~ T ' (3) 
where 7 is the Lorentz factor of the particle, 7 = (1 — t^/c 2 )" 1 / 2 . In the target frame, this 
equation applies, as long as 7 is taken as the boost to go from the frame of one nucleus to 
the other (7 = 2^ ollider - 1). 

For a grazing collision, where the two nuclei barely touch, we can take b min = 2Ra, and 
the maximum photon energy is r yhv/2R^ (Ra = nucleus radius). The maximum photon 
energy is about h/ (2RAAm p c) of the ion energy. Here, Am p is the ion mass. For heavy-ions, 
Ra ~ 7 fm so oj max is about 0.03/A of the ion energy. For protons, Ra is not well denned, 
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TABLE I: Some ion species, maximum energy and luminosity for several accelerators [17PJ. Also 
shown are the maximum effective 7p and 77 energies. For proton beams, the maximum effective 
photon energy is taken to be 10% of the proton energy, although there is some flux at higher 
energies. The CERN SPS is a fixed target accelerator; the effective luminosity depends on the target 
thickness. Not mentioned here are lower-energy accelerators, where photon exchange processes have 
also been studied. 

Accelerator Ions Max. Energy (CM) Luminosity Max. 7A Max. 77 
per Nucleon pair (cm _2 s _1 Energy Energy 



CERN SPS Pb+Pb 


17 


GeV 




3.1 


GeV 


0.8 


GeV 


RHIC Au+Au 


200 


GeV 


4 x 10 26 


24 


GeV 


6.0 


GeV 


RHIC p+p 


500 


GeV 


6 x 10 30 


79 


GeV 


50 


GeV 


LHC Pb+Pb 


5.6 


TeV 


10 27 


705 


GeV 


178 


GeV 


LHC p+p 


14 


TeV 


10 34 


3.1 


TeV 


1.4 


TeV 



but taking uj max to be 10% of the proton energy is a reasonable rule-of-thumb. RHIC (see 
Table 1) can reach photon-gold center of mass energies up to about 30 GeV per nucleon, 
and photon-proton center of mass energies up to 300 GeV. These energies are slightly higher 
than are available at fixed target accelerators and at HERA respectively. At the LHC, the 
corresponding energies are 1 TeV and 10 TeV respectively, more than an order of magnitude 
higher than is available elsewhere. 

The equivalent (or virtual) photon flux per unit area (the relation between n{u) and 



N(u,b) is n{uj) = f N(u,b)d 2 b) is |]| 



17] 



where x = ub/ 'yfihc, Z is the ion charge, a = 1/137, [3c is the particle velocity and K and K\ 
are modified Bessel functions. The first term (Ki(x) 2 ) gives the flux of photons transversely 
polarized to the ion direction and the second is the flux for longitudinally polarized photons. 
The transverse polarization dominates for ultra-relativistic particles (7 ^> 1). The photon 
flux is exponentially suppressed when u > ^f3hc/b, justifying the estimates in the beginning 
of this section. 
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These photons are almost real, with virtuality — q 2 < (H/ ' Ra) 2 '■ Except for the production 
of e + e~ pairs, the photons can usually be treated as real photons. 

The usable photon flux depends on the geometry. Most UPC reactions lead to final 
states with a handful of particles. These final states will be overwhelmed by any hadronic 
interactions between the fast moving ion and the target. Thus, the useful photon flux is that 
for which the ions do not overlap, i.e. when the impact parameter b = | t»i — b>2 1 is greater 
than twice the nuclear radius (2Ra) (see Figure ^c)). Usually, we can take R A = 1.2A 1 ' 3 
fm, where A is the atomic number. The b > 2Ra requirement treats the nuclei as hard 
spheres; it is accurate for heavy nuclei, but less appropriate for lighter ions. 

The photons can interact with a target nucleus in a one-photon process, (when bi < Ra) 
or with its electromagnetic field in a two-photon process when bi > Ra and b 2 > Ra- In 
a photonuclear (one-photon) interaction, the usable photon flux is obtained by integrating 
Equation |U over b > b m i n = 2Ra'- 



2Z 2 a 



n[uj) 



2 



^o(0^i(0-|-(^i 2 (0-^o 2 (0) 



(5) 



7Tf3 

where £ = u)b min / '7 (3hc = 2uRA/ r fPhc. 

For two-photon exchange processes, the equivalent photon numbers, in Equation El must 
account for the electric field orientation of the photon fluxes with respect to each ion (see 
Figure ^b)), obeying the ion non-overlap criteria 61, 6 2 > Ra Ql- The field orientation is 
not included in Equation |S| For instance, owing to symmetry properties, J n = + (scalar) 
particles originate from confi gur ations such that Ex || E 2 , whereas 0~ (pseudo-scalar) parti- 
cles originate from E x _L E 2 jisl l"!^. If one uses Equation for n(ui) and n(u; 2 ), the total 
photoproduction cross section obtained from Equation |2] is higher than in a more detailed 
calculation, and the difference increases with increasing particle masses Q|. Even more 
detailed calculations can be done by replacing the sharp-cutoff, b\, 6 2 > Ra, criterion with 
integrals over hi and b 2 which are weighted by the hadronic non-interaction probability. 
Asymmetric collisions (especially pA and dA) are also of interest; the higher-Z nucleus is 
likely to be the photon emitter, so the photon direction is known. 

Low energy processes, e.g. nuclear excitation are also sensitive to the electromagnetic 
multipolarity involved. Equations |U and El are only appropriate for electric dipole (El) 
excitations. Equations for higher multipolarities are described in Reference 0. 

For protons, the hard sphere approximation is inadequate. Instead, the proton size is 
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included by the use of a form factor. With a dipole form factor, the flux is 



n(uj) 



2ttz 
where 



[l + (1 _ 2)2] ( lnx _ii + i__L + _L) (6) 

accounts for the proton structure and z = W 2 /s, with W the 7p center of mass energy, 
and s the squared ion-ion center of mass energy per-nucleon. Here, Q m in is the minimum 
momentum transfer possible in the reaction. For proton-proton collisions, the form factor 
has an effect similar to imposing a requirement b m i n = 0.7 fm 21 1. 



For protons and light nuclei, the weak electromagnetic interactions introduce another 
complication. The momentum transfer due to elastic scattering, Ap = 2rjh/b, with rj = 
Z\Ziol <d 1 is small enough that the impact parameter is not a well defined observable 
(because ApAb ~ H, leading to Ab > b for rj < 1) 22]. This does not affect the total photon 
flux. However, it might affect the component of the photon flux that is unaccompanied 
by hadronic interactions. The uncertainty may also affect the probabilities for multiple 
interactions, discussed in Section 5. 

Equation El is valid when the proton remains intact. When photon emission with proton 
excitation, such as to the A resonance, is included, then the flux increases about 30% [23]. 



At very hig. 
important 



iphoton energy (z — » 1), the magnetic form factor of the proton can also become 



B. Experimental Characterization 

Ultra-peripheral collisions look very different from the more conventional hadronic in- 
teractions. The final state multiplicity is much smaller, and, usually the events are fully 
reconstructed. Because the photon p? are small (p? ~ lu/^c), the final state px will also be 
small. Photonuclear interactions that involve coherent scattering from the target nucleus 
(such as vector meson production) also have a very small px'- Pt < h/ Ra- This gives the 
events a distinctive experimental signature, greatly simplifying detection 25| . 

UPCs are studied at a variety of accelerators. The characteristics of some relevant accel- 
erators are given in Table |U Each accelerator can accelerate many different species; Table H] 
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gives only a few candidates. The CERN Super Proton Synchrotron (SPS) has produced 
results on lead-to-bismuth transmutation and e + e~ pair production in ion-ion collisions. 

Although RHIC only began taking data in 2000, it has already released UPC results on 
p° photoproduction and on e + e" pair production. RHIC has enough energy and luminosity 
to photoproduce a wide variety of mesons, including the J/if>. However, because it is a 
collider, detection of very low pr particles is difficult, complicating the study of e + e~ pairs 
and other atomic phenomena. 

Although it is exclusively a pp collider, the Fermilab Tevatron is an interesting place 
to study UPCs. Antihydrogen was produced there using the process 77 — > e + e~, with 
the positron bound to an antiproton |a|. Photoproduction of the J/ip 21 1 may have been 
observed by the CDF collaboration 26j] . 

The Large Hadron Collider (LHC), scheduled to begin operation in 2007, will search for 
physics beyond the standard model. A UPC program at the LHC can contribute to this 
search. Especially for pp collisions, where 77 and 7p energies up to about 10% of the beam 
energy are accessible, UPCs may an attractive place to search for new physics. With ion 
beams, the photon energies are lower, but W ± , Z, and heavy quark physics may be studied. 



II. LOW ENERGY PHOTONUCLEAR INTERACTIONS 



Relativistic Coulomb excitation (RCE) is now a popular tool to investigate the intrinsic 



nuclear dynamics and structure of the co 
involving radioactive nuclear beams J27, 



liding nuclei. It is especially important in reactions 



28 



29 



30 



31 



32, 



and has been used for many 
decades in low energy nuclear collisions to study nuclear structure 34]. However, nuclear- 
induced processes may also contribute to the reactions being studied. 

RCE may involve single or multiple photon-exchange between the projectile and the tar- 
get. In the first case, perturbation theory directly relates the data to the matrix elements 
of electromagnetic transitions. These matrix elements are clean probes of the nuclear struc- 



ture, and RCE can be used to study short- 

nr 

real photons or electron scattering [28J, |31 , 



ived unstable nuclei that cannot be probed with 



33] 



Radiative capture processes (b + c — ► a + 7) play a major role in astrophysical sites, e.g., 

Some reactions of interest for astrophysics, e.g., 7 Be (p, 7) 8 B, 



m a pre- supernova 



can be studied via the inverse photo-dissociation reaction 8 B (7,p) 7 Be |3j| using relativistic 
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FIG. 3: (a) S-factors (Si?) forjdie 7 Be(p,7) 8 B reaction. The GSI data was obtained using the 
Coulomb dissociation method 



■*'5<S]. The oilier dala arc Iron) direct capture measurements 
(b) Cross sections for the excitation of the GDR (1-phonon) and the DGDR (2-phonon) in 208 Pb 
projectiles incident on different targets. The dashed curves are theoretical calculations. 

Coulomb collisions. The Coulomb breakup reaction a + A — > b + c + A is useful to obtain 
the corresponding 7 induced cross section 7 + a — > b + c. Using detailed balance, this 
cross section can be related to the radiative capture cross section b + c — ► a + 7, of 



astrophysical interest [37[ . The radiative capture cross sections are often expressed in terms 
of the astrophysical S-factor: S(E) = <j{E) exp(—2iTZbZ c e 2 /hvb c )/ E, where E = E re \ = 
m bc vl c /2 is the relative energy between b and c. In this equation t> fec is the relative velocity 
and rribc = TnbTn c / (nib + Tn c ) is the reduced mass of b + c. Because the Coulomb penetration 
factor is explicitly factored out, the S-factor is a much natter function of E than u(E) 
allowing a better extrapolation of the measurements. 

As an example, Figure Ufa) shows the result of an experiment performed at the GSI 
laboratory, in Darmstadt, Germany [3| for the Coulomb dissociation of 8 B. Data on the 
reaction 7 Be + p — > 7 + 8 B is important for understanding the structure of our sun. The 
decay of 8 B is responsible for the high energy neutrinos observed by earth-bound detectors. 
The measured S-factor (S17, l=proton, 7= 7 Be) is shown in Figure Efa) as solid circles. The 
solid curve is a fit using a theoretical model for S 17(E). Some of the data shown in the figure 
are from direct capture experiments 
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Other b(c, 7)0 radiative capture reactions are planned, or have already been studied with 

n 

the Coulomb dissociation method 40J. These processes may occur in the sun, supernovae, 
or during the Big Bang. Most of these reactions cannot be directly studied because the 
Coulomb barrier leads to very small values for the cross section, beyond the reach of present 
experimental techniques. 



Giant dipole resonances (GDR) occur in nuclei at energies of 10-20 MeV. Their gross 
properties are well described in terms of an out-of-phase collective motion (oscillation) of 



protons against neutrons in a nucleus 



41 



, . If this oscillation is harmonic, with excitation 
energy ftw, then higher excitation modes, with energies equal to Nhw also occur. These 
modes are interpreted as double, triple, giant dipole resonances. The double giant dipole 



resonances (DGDR) are thus two giant dipole vibrations superimposed in one nucleus, with 



about twice the energy of the GDR 



nt dipole vi br 



In the harmonic model, the RCE cross 



section for all multiphonon states can be calculated exactly jj] 

A series of experiments at the GSI laboratory obtained the energy spectra, cross sections 



44, 



47 



and angular distribution of fragments following the decay of the DGDR [42 
Q Q]- The experimental cross sections are about 30% bigger than the theoretical ones. 
This is shown in Figure 1Mb) where the dashed lines are the result of theoretical calculations 



5 it l52t l53t l54j . These experiments are promising for the studies of the nuclear 
response in very collective states. The giant resonances are still poorly understood, even 
with the best current microscopic approaches. The study of the width of the DGDRs will 
be helpful to improve this scenario 30]. 

In heavy ion colliders the mutual Coulomb excitation of the ions (leading to their simul- 
taneous fragmentation) is a useful tool for beam monitoring 55]. A recent measurement at 
RHIC Q, 

using the Zero Degree Calorimeters to measure the neutron decay of the reaction 
products, has proved the feasibility of the method. 

DGDRs constitute only 10% of the total fragmentation cross section induced by Coulomb 
excitation in UPC. The dominant contribution is the excitation of a single GDR, which then 
decays mostly by neutron emission. This is also a major source of beam losses in relativistic 
heavy ion colliders 57J , and an important fragmentation mode of relativistic nuclei in cosmic 
rays. 

Another useful reaction is deuteron photodissociation in d+A collisions - a photon from 



11 



a heavy-ion photodissociates a deuteron 58]. The reaction has a large cross section, 1.38 
b for d+Au at RHIC, and 2.49 b for d+Pb at the LHC |59j, and has been used as a 
'standard candle' for luminosity monitoring 
are technically simpler than p+A collisions. 



6Cj. d+A collisions are studied because they 



III. PHOTOPRODUCTION AT HADRON COLLIDERS 

The main interest in photoproduction at hadron colliders is derived from the possibility it 
offers of a direct determination of the gluon distribution in nucleons and nuclei. Examples of 
interactions in which the gluon distribution can be probed are exclusive production of heavy 
vector mesons, photoproduction of heavy quark- ant i- quark pairs, and photoproduction of 
jets. These gluon distributions are not directly accessible in deep inelastic scattering, because 
the gluons carry neither electrical nor weak charge. 

Measuring the nuclear shadowing using heavy-ion beams is particularly interesting. The 
nuclear gluon density first approximation, be written as the nucleon gluon distri- 

bution, g(x,Q 2 ), multiplied by the number of nucleons (A): 

G A (x,Q 2 )=Ag(x,Q 2 ). (8) 

Here, x is the fraction of the projectile momentum carried by the gluon, and Q 2 is the 
4-momentum transfer squared. 

Results from deep inelastic scattering of electrons on nuclear targets have, however, 
showed deviations from such a simple scaling for the structure function, F 2 (x 2 , Q 2 ). Depend- 
ing on x and Q 2 , suppression (shadowing) of up to ~30% and enhancement ( ant i- shadowing) 
of up to ~10% have been observed. The effects of shadowing on G A (x, Q 2 ) are hard to de- 
termine directly. The current best estimates of the modification to the gluon distribution in 
nuclei are obtained from the Q 2 -evolution of F 2 (x 2 , Q 2 ) 61, and from studies of diffrac- 
tive interactions Photoproduction at heavy-ion colliders may provide a more direct 

measurement of G A (x,Q 2 ). 

The particle production in photon-hadron or photon-nucleus interactions can be exclu- 
sive, when the protons or nuclei remain in their ground state or are only internally excited, 
or inclusive, when at least one of the nucleons or nuclei breaks up. Exclusive production 
will be discussed first. When the momentum transfer is small compared with its inverse 
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nucleon/nuclear size, Q ~ he/ R, the fields couple coherently to the entire target. The kine- 
matics of coherent, exclusive interactions is very similar to that of two-photon interactions, 
which will be discussed in section 4. 



A. Exclusive Particle Production 

The dominant coherent interaction leading to the production of a hadronic final state is 
the exclusive production of vector mesons, 

A + A^A + A + V. (9) 

In these reactions a photon from the electromagnetic field of one of the projectiles interacts 
coherently with the nuclear field of the other (target) producing the vector meson. 

Exclusive vector meson photoproduction on proton and nuclear targets has been studied 
since the mid 1960's using photon beams and since 1992 at the HERA electron-proton 
accelerator j^. The first results from a heavy ion collider on exclusive p° production 
Au + Au — > Au + Au + p°) were recently published by the STAR collaboration at RHIC 

The total vector meson cross section in p+p or A+A interactions can be calculated 
from Equation By differentiating and changing variable from uj to y, the rapidity of the 
produced vector meson, one obtains 

da(A + A^A + A + V) 

where the photon energy, uj, is related to y through u = (M v c 2 /2) exp(y) and My is the 
mass of the vector meson. If the photon flux is known, the differential cross section, da/dy, is 
thus a direct measure of the vector meson photoproduction cross section for a given photon 
energy. 

The bulk of the photon-hadron cross section can be explained by the photon first fluctu- 
ating to a qq pair, which interacts with the target through the strong nuclear force. Since 
the photon has quantum numbers J PC = 1 , it preferentially fluctuates to a vector meson. 
The lifetime of the fluctuation is determined by the uncertainty principle. For a photon of 
virtuality Q fluctuating to a state of mass My the lifetime is of the order of 

^ ^ /-.-.X 

^ ~ = W — — . (11) 

y/M y C 4 + Q 2 C 2 MyC 2 
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The last approximation holds at hadron colliders because of the low virtuality of the photons. 
The photon wave function is written as a Fock decomposition 6(|: 

\-i >=C ba re\lbare> +C p \p> +C U \U0> +C 4 ,\(j)> +...+C q \qq> . (12) 

Here C bare ~ 1 and Cy ~ ^/a em (V = p, 00, <f), ■ ■ ■ ). The coefficients Cy are related to the 
photon- vector meson coupling, fy, through 



Cy = (13) 
JV 

The numerical values of the couplings fy are usually determined from the vector meson 
leptonic decay widths, T{V — > e + e~). 

According to the Generalized Vector Meson Dominance Model (GVMD), the scattering 
amplitude for the process 7 + A — > B is the sum over the corresponding vector meson 
scattering amplitudes, 

A 1+A ^ B {s, t)=J2 C v A v+A ^ B {s, t) . (14) 



v 



For "e 



astic" scattering, 7 + A — > V + A, the cross-terms, i.e. V + A — > V + A, are usually 
and are often neglected. The cross section is then (t is the momentum transfer 



from the target nucleus squared and da/dt = \A 



2\ 



da^ + A^V + A) ^ 2 da{V + A^V + A) 

It = ° v * ' (15) 



The momentum transfer of the elastic scattering is determined by a hadronic form factor, 
F(t), 

\F(t)\*. (16) 



da da 
dt dt 



2 



For proton targets, the form factor is well represented by an exponential function, |-F(t)|" 
exp(— b\t\) with a slope b w 10 GeV -2 c 2 for the light vector mesons (p, to) and b ~ 4 GeV -2 c 2 
for the J /ty. The form factor for nuclear targets is peaked at much smaller momentum 
transfers because of the larger size of the target. 

The form factor reflects the size and shape of the target. It can in principle be calculated 
if the spatial distribution is known. The dynamical information is contained in the forward 
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TABLE II: Cross sections for exclusive vector meson production in Au+Au and Pb+Pb interactions 
at RHIC and the LHC, respectively j^ . 



Meson Au+Au, RHIC Pb+Pb, LHC 





a [mb] 


a [mb] 


p° 


590 


5200 


OJ 


59 


490 




39 


460 




0.29 


32 



scattering amplitude, da/dt(t = 0). The optical theorem relates this to the total vector 
meson cross section, a to tiyA): 



da 
~dt 



t=0 



lQirh 2 



(17) 



Here, rj is the ratio of the real to the imaginary part of the scattering amplitude. 

In Reference j68j, data on vector meson photoproduction with proton targets were used 
to extract the total vector meson nucleon cross section, a tot (VN). This was then used to 
calculate the total vector meson nucleus cross section, a to t(yA), from the nuclear geometry. 
This gave the vector meson production cross sections for heavy ion interactions at RHIC 



and the LHC shown in Table 
may reduce the cross-section 



For heavier vector mesons, like the J/ip, gluon shadowing 



In the Glauber Model jz(|, the elastic scattering amplitude is given by the two-dimensional 
Fourier transform of the nuclear profile function, T(b): 

2 



da(^ + A^V + A) _ 7T 

It = ¥ 



J e iPT-b/fc r ( b ) d 2 h 



(18) 



T(b) is a function of the distribution of matter inside the nucleus, p(b, z), and the vector 
meson-nucleon forward scattering amplitude, fvN (which can be related to the total vector 
meson-nucleon cross section through Equation fTTj) : 



r(b) 



exp 



U3 



p(b,z')f VN (0)dz' 



(19) 
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This approach only works for high photon energies, when cry [3 At > Ra so the interaction 
is longitudinally coherent over the entire nucleus. At lower u, the loss of coherence reduces 

n n 

the cross section. The Glauber model is discussed in References |7Jj and |72j . 

A Glauber model calculation of the coherent p° production cross section in Au+Au 

n 

collisions at RHIC gave a total cross section of 934 mb [73j. This is about 50% higher 



than the result in 



681 ] (cf. Table ITT]) . The main reason for the difference is that in 68] the 



total vector meson-nucleus cross section was calculated assuming that a to t(pA) ~ <j ine i(pA). 
The calculation in 



|73j furthermore includes the contribution from off-diagonal elements 
corresponding to p' + Au — > p + Au scattering, as well as a non-zero real part of the forward 
scattering amplitudes (77 in Equation [T7J). For a discussion of the p'+A — > p+A contribution, 
see also [67]. 

The measured p° production cross section at RHIC is cr(Au + Au — ► Au + Au + p°) = 



460 ± 220 ± 110 mb at 130 A GeV l25j. This can be compared with the Glauber Model 

n 

calculations, which give a = 490 mb at this energy [74]. The corresponding number from 
the method used in [6.8] is o = 350 mb. 

The rapidity distribution for coherent p° production measured by the STAR collaboration 
in Au+Au interactions at 200 A GeV is shown in Ufa) [75] . This is for p° production in 
coincidence with mutual Coulomb breakup of the beam nuclei (cf. Section 5). The rapidity 
distribution and cross section are in excellent agreement with the distribution obtained from 
the Monte Carlo m odel based on the captions in 0, coated for the e„ta, 
acceptance. The reconstructed invariant 7t + tt~ mass is shown in Figure |U (b). The shape is 



ativistic Breit-Wigner function and a Soding interference 



well described by the sum of a re 
term for direct 7r + 7r~ production 



The STAR collaboration has very recently presented the first preliminary data on p° 
production in deuteron-gold interactions (the gold nucleus acts as photon emitter) |[[8| and 
on coherent production of 7r + 7r~7r + 7r~ [79]; the latter may be attributed to p* photoproduc- 
tion. PHENIX has shown indications of coherent J/ty and e + e~-pair production in Au+Au 
interactions at RHIC [goI . Isoj ] . 

The forward scattering amplitude for heavy vector mesons has been calculated from two- 
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FIG. 4: Rapidity (a) and invariant mass (b) distributions for coherent p° production in Au+Au 
interactions accompanied by mutual Coulomb breakup at y/s = 200 A GeV, by the STAR collab- 
oration. The dashed curves in b) corresponds to a relativistic Breit-Wigner function and a Soding 
interference term; the solid curve is the sum of the two. The dash-dotted curve describes the 

n 

background from incoherent interactions 175!] . 

gluon exchange in QCD. To leading-order 
da{^p — > Vp) 



dt SeSf^Nl-.^/")] 2 - (20) 

Here, x is the fraction of the proton or nucleon momentum carried by the gluons and 
the gluon distribution, g(x,Q 2 ), is evaluated at a momentum transfer Q 2 = (My/2) 2 . This 
approach has been developed further by including relativistic wave functions and off-diagonal 



parton distributions |8J 



831 ] . The result is a total vector meson nucleon cross section which 



grows rapidly with increasing photon-proton center-of-mass energy, W lp . For T production, 
cr oc is expected. 

The dependence of dojdt on [g(a;)] 2 makes exclusive vector meson production a very 
sensitive probe of the proton and nuclear gluon distributions. An T-meson produced at 
mid-rapidity at the LHC would come from gluons with x ~ 7 ■ 10~ 4 and x ~ 2 • 10~ 3 in p+p 
and Pb+Pb interactions, respectively. 

Figure shows the predicted da/dy for heavy vector mesons in nucleus-nucleus and 
proton-proton collisions. The calculations are based on parameterizations of the photon- 
oroton cross sections derived from measurements at HERA and from QCD based models 
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68). 
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RHIC Au+Au -> Au+Au+J/ 1 ? 



LHC Pb+Pb -» Pb+Pb + T 




FIG. 5: Rapidity distributions for exclusive J/ip and T production in nucleus-nucleus and proton- 
proton collisions. Adapted from and j^ . 



The photoproduction cross sections rise rapidly with energy. With lead at the LHC, the 
p° cross section is comparable to the hadronic cross section. With calcium beams at full 
luminosity, the LHC will produce about 230,000 p°, 15,000 <fi and 800 J /ip per second. These 
rates are comparable to dedicated e + e" colliders, qualifying the LHC as a meson factory! 

B. Interference in Exclusive Vector Meson Production 



When a single vector meson is produced through a coherent photonuclear interaction in 
a nucleus- nucleus or proton- (anti-)proton collision, it is in general not possible to determine 
which projectile acted as target and which was the photon-emitter. The two possibilities are 
indistinguishable and under certain conditions they will interfere quantum mechanically [84] . 
Because of this interference, it is incorrect to add the cross sections for the two possibilities. 
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The cross section is given by adding the corresponding amplitudes A\ and A 2 



hda 



A 1 ±A 2 \ 2 d 2 h. (21) 



dydpr j b>2R 

The interference is maximal at mid-rapidity, where symmetry requires that A\ — A 2 . For 
ion-ion and proton-proton collisions, the interference is destructive because of the negative 
parity of the vector meson; exchanging the position of the two nuclei or the two protons 
is equivalent to a reflection of the spatial coordinates, i.e. a parity transformation. For pp 
collisions, as at the Fermilab Tevatron, exchanging the proton and antiproton involves a 
charge-parity (CP) transformation. Since CP is positive for vector mesons, the interference 



is constructive at pp colliders 



2l|. 



The amplitudes A\ and A 2 depend on the photon flux (e.g. rapidity) and on the pho- 
tonuclear cross sections. Their px dependence comes from the convolution of the photon px 
spectrum and the pt from the photon-nucleus scattering. The former is given by the equiv- 
alent photon spectrum 0,0], and the latter comes from the form factor of the target. 

If the outgoing vector meson is treated as a plane wave (appropriate for a distant ob- 
server), at mid-rapidity, A\ = A 2 and the square of the sum of the amplitudes is 

\A 1 ±A 2 \ 2 = 2J^(l±cos(^-)\ . (22) 

For very low momenta, px <C cos(p • b/h) « 1 — (p • h/h) 2 /2 and, as px — > 0, 

the interference is complete; emission disappears in ion-ion collisions, but doubles for pp 
colliders. Interfence is significant for px < 20 MeV/c for the p° at RHIC |8J], and px < 250 
MeV/c for the J/ip at the Tevatron [2l| . When b ^> h/px, the cosine term oscillates rapidly 
as b varies, and the interference disappears. In this regime, the cross section reduces to the 
sum of cross sections for the two photon directions. 

Away from mid-rapidity, \Ai\ ^ \A 2 \ because the photon energies for the two possibilities 
are different: o; 1>2 = (M v c 2 /2) exp(±y). Both the photon flux and photonuclear cross 
sections will be different. The interference will thus be reduced. A\ and A 2 could also 
have slightly different phases, adding a phase factor S to the cosine term in Equation l2*2l 
However, 5 = in the standard Pomeron models Q|, and a significant phase difference 
seems unlikely. 

The interference in exclusive vector meson production is of particular interest because 
the vector mesons have very short lifetimes compared to the typical impact parameters. The 
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FIG. 6: Efficiency corrected t± spectrum for p° production in Au+Au collisions at RHIC with 
Coulomb breakup. The points are the data, and the solid curve is a fit to Equation 1231 Interference 



causes the dip at low t < 0.001 (GeV/c 



,2 



median impact parameter for p° production in Au+Au collisions at RHIC, for example, is 
46 fm, much larger than the lifetime of the p° (r = 1.3 fm/c). The vector meson cannot, 
on the other hand, be produced more than ~1 fm away from one of the two nuclei because 
of the short range of the nuclear force. Observing the expected interference pattern would 
thus prove that the wave function of the vector meson is preserved long after it has decayed. 
This is an example of the Einstein- Podolsky- Rosen paradox js^J. 

Preliminary data from the STAR collaboration on p° production in Au+Au at y/s nn = 
200 GeV seem to confirm the presence of interference 86j . The measured t spectrum is shown 
in Figure |3 The data are for interactions where both gold nuclei Coulomb dissociate. The 
coincident Coulomb dissociation selects events with smaller impact parameters compared 



to exclusive production 
parameter is only 18 fm 



cf. Section 5); with Coulomb dissociation, the median impact 



The data are fit to a function, 

dN 

— = a exp(-6|t|) [1 + c(R(t) - 1)] , (23) 

with three parameters. These correspond to a normalization constant (a), the width of 
the nuclear form factor (b « R\/h 2 ), and a parameter to quantify the magnitude of the 
interference (c). The function R(t) is the ratio of Monte Carlo da/dt calculated with and 
without interference. This functional form separates the interference from the nuclear form 
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factor. No interference would correspond to c = 0, while complete interference according to 
the calculations above would correspond to c = 1. A fit to the data finds c = 1.01 ± 0.08 
(0.1 < y < 0.5) and c = 0.78 ± 0.13 (0.5 < y < 1.0) for the two ranges in rapidity. 



C. Inclusive Photoproduction 

The high photon flux at hadron colliders and the large total photon-hadron cross sections 
lead to high rates for other photonuclear interactions. In Au+Au collisions at RHIC, the 
total photonuclear cross section for photon-nucleon center-of-mass energies above 4 GeV is 
about 2 barns, or nearly 1/3 of the total hadronic Au+Au cross section. The majority of 
these interactions are resolved interactions, i.e. they are preceded by a fluctuation of the 
photon to a qq state. They therefore resemble inelastic hadron- nucleon/nucleus collisions. 
Because the photon energies are much lower than the beam energies, the kinematics is similar 
to that of fixed target interactions. 

Despite the large cross section, particle production in resolved photon- nucleon/nucleus 
interactions at hadron colliders has so far attracted relatively little interest. See, however, 
Reference Understanding the kinematics of these interactions is nevertheless essential, 
because they are a significant background to other ultra-peripheral processes, particularly 
at the trigger level 88]. 

Considerably more interest has been devoted to direct photon interactions, in particular 



the production of heavy gg-pairs 



Recently, the cross section for photoproduc- 



tion of heavy quark pairs in pp collisions has been calculated |92|. In these interactions, a 
photon interacts with a parton in the target and the partonic cross section can be calculated 
from QCD. 

The leading order contribution to the photoproduction of a gg-pair corresponds to photon- 
gluon fusion, as is illustrated in the Feynman-diagram in Figure H(d). The cross section for 



the partonic sub-process is 



93, 



94] 



(T lg ^ m {W 1 



ne 2 a em a s (Q 2 )h 2 c 2 



'■yg^qq V" 79/ ' ti/2 

79 



(3 - /5 4 ) ln(^|) - 2(3(2 - (3 2 ) 



(24) 



Here, m q and e q are the quark mass and electric charge, respectively, (3 = (1 — 4m^c 4 /W 2 g ) 
and W ig is the photon-gluon center-of-mass energy. If the gluon carries a fraction x of the 
nucleon momentum, then W 2 g = 2uxy/s. The strong coupling constant, a s (Q 2 ), is evaluated 
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FIG. 7: Rapidity distributions of bottom (left) and top (right) quarks produced in photonuclear 
collisions at the LHC. The solid and dashed curves are for Pb+Pb and O+O interactions, respec- 
tively. Here, the photon is emitted by the nucleus with positive rapidity; the complete cross section 
is the sum of this curve plus it's mirror image. Shadowing is included using the parameterizations 
of Reference |62j ] . Adapted from References 

Q and 0- 

to one loop at scale Q 2 = m 2 c 2 + pf,, where px is the quark transverse momentum. 

The total photoproduction cross section a (A [7] .A — > AqqX) is obtained by convoluting 
the partonic cross section with the equivalent photon flux, n(u), and the nuclear /nucleon 
gluon distribution, G A (x2,Q 2 ), i.e. 

a(A[j]A -> AqqX) = J J G A {x, Q 2 ) a ig {W ig ) ®(W ig - 2m q c 2 ) dudx . (25) 

This equation is the equivalent of Equation |21 for two-photon interactions with the photon 
flux from one nucleus replaced by the gluon distribution, G A (x, Q 2 ). The final state qq 
rapidity depends on the photon energy and the gluon x. The rapidity distributions of 
bottom and top quarks produced in Pb+Pb and O+O collisions at the LHC are shown in 
Figure [7| The kinematics are discussed in more detail elsewhere 91 



The qq production cross section is peaked near threshold, W ig w Am 2 . Mid-rapidity 
production of cc- and 66-pairs therefore mainly probes x-values of x ~ 1 • 10~ 3 (cc) and 
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TABLE III: Cross sections for qq photoproduction through direct photon- gluon fusion in heavy ion 
interactions. The numbers in column 4 and 5 include nuclear gluon shadowing from References 
jo^ and 97], respectively. 



qq cross sections in heavy ion collisions 



Colliding system 


flavor 


a [mb] 


a [mb] 


a [mb] 






No shadowing 


EKS98 


FGS 


RHIC Au+Au 


cc 


15.8 


17.4 


17.6 


RHIC Au+Au 


bb 


2.9 -1CT 3 


3.0 d0~ 3 


3.0 -10~ 3 


RHIC Si+Si 


cc 


0.196 


0.203 


0.20 


RHIC Si+Si 


bb 


1.07 -10~ 4 


1.13 -10~ 4 


1.14 -10~ 4 


LHC Pb+Pb 


cc 


1250 


1050 


850 


LHC Pb+Pb 


bb 


4.9 


4.7 


4.4 


LHC Ar+Ar 


cc 


16.3 


14.3 


12.3 


LHC Ar+Ar 


bb 


0.073 


0.070 


0.066 



x ~ 3 • 10 3 (bb) in heavy-ion collisions at the LHC. The corresponding numbers at RHIC 
are x ~ 2 ■ 10~ 2 and x ~ 1 • 10 _1 . (For a gg-pair with invariant mass W 7g and pair-rapidity 

y, x = (W 7fl A/i)e*.) 

The total cross sections for cc and bb production in various systems at RHIC and the 
LHC are listed in Table IIHI The calculations without shadowing are compared with two 
calculations that include nuclear modifications. As expected, shadowing has the largest 
effect on the production of lighter quarks (cc) using heavy nuclei. The cross section for 
producing a cc-pair in Pb+Pb interactions at the LHC is of the order of 1 b, nearly 1/6 of 
the total hadronic cross section. 



Quark pairs can also be produced in anomalous interactions, where a parton from the 
resolved photon interacts with a parton in the target, or in two-photon interactions. The 
cross sections from anomalous interactions are small compared with the direct production 
cross sections j^l|. The anomalous cross sections are 1-20% of the direct cross sections, 
depending on quark flavor and collision energy. The two-photon contribution is usually less 
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than 1% of the anomalous cross section. 

In addition to probing the nuclear gluon distribution, the photon-gluon fusion reactions 
are of interest as a means to determine the electric charge of the top quark. The top quark 
measurements are all consistent with a standard model top with electric charge q t = +2/3. 
However, since the correlation between the decay products W + «-> b or W~ «-> b are never 



measured, q t is unconstrained 



961 ] . The data can still accommodate an exotic particle t' 



with q t = —4/3, which decays via t' — > W~ + b. Since the cross section for 7 + g— > q + q is 
proportional to qf , q t = -4/3 would quadruple the cross section. 

There have so far been no experimental measurements of heavy quark pair photoproduc- 
tion in heavy ion or pp collisions. Some of the experimental techniques that could be used 
to separate this signal from hadronic backgrounds are discussed in Section 6. 



D. Dijets, Compton scattering, Vector Boson Production and other processes 

When the final state quarks in the process shown in Figure Q (d) have high p?, the final 
state is two roughly back-to-back jets 98[. The cross section is calculated in a manner 
similar to that of Equation |2U except that light quarks are included. Jets may also be 
produced to leading order through the process 7 + q — > g + q. The cross section for dijet 
production is sensitive to the gluon distribution in the target; because of the simplicity of 
the reaction, there are fewer systematic uncertainties than in other processes, like vector 
meson production. However, it may be more difficult to isolate dijet photoproduction from 
backgrounds, such as hadronic production of dijets and diffractive hadronic production. 

Since the gluon-contributing nucleus does not stay intact, the experimental signature 
for this process is two jets, accompanied by a single rapidity gap between the jets and 
the photon-emitting nucleus. The two jets may have very different rapidities and it may 
be difficult to reconstruct the entire event. Calculations have considered the case where a 
single jet is detected, with \rj\ < 1 (9^. Without shadowing, the rate to photoproduce jets 
with energies above 21 GeV in lead-lead collisions at the LHC is 0.015 Hz. In a 10 6 s run, 
jets up to 80 GeV should be detectable. 

A closely related process is the production of a photon + jet final state; this is essentially 
Compton scattering. The rates for this process are about two orders of magnitude below 
that for dijet production 
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The strong Coulomb fields may also dissociate hadrons into jets. For example, a proton 
may fragment into three quarks, leading to reactions such as 7p — > 3 jets; this would be a 



distinctive signature in pA collisions 
to reactions such as 77 — > 2 jets (3, 



99]. One photon can also dissociate another, leading 
IOC 

W ± and Z° can be photoproduced in ultra-peripheral collisions. A Z° can be produced 
when a photon fluctuates to a qq pair, scatters from a target nucleus an d em erges as a Z°. In 
the high-energy limit, the cross section for 7p — ► Z°p is about 0.01 pb 10l| . Unfortunately, 
even with a coherent photon beam, the cross section seems too low to be observable. 

pairs can be produced directly from a photon fluctuation, 7 — > W + W~ . One of 
the Ws can interact with the target nucleus, leading to a hadronic jet plus a real W. 
Unfortunately, this process has not been studied in detail. 



IV. TWO-PHOTON PROCESSES 



A. Production of free- and bound-pairs 



Between 1933 and 1937, Furry, Carlson, Landau, Lifshitz, Bhabha, Racah, Nishina, 
Tomonaga, and several others perfor med calcu l ations of e + e~ production in relativistic col- 



102 



103 



104, 



105 



106]. The purpose was to test the 



lisions of fast particles (cosmic rays) 
newly born Dirac theory for th e pos itron. Starting with the Dirac equation for the electron 



and its antiparticle they found 
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a = ct [L 3 - 2.198L 2 + 3.821L - 1.632] 



(26) 



27tt 

where <To = (ZiZ2a 2 h/m e c) 2 , L = 111(7172), and 7^ is the Lorentz factor of ion % in the 
laboratory system. The first term of this equation can be simply obtained from Equation 
1 and the cross sections for 77 pair-production. The production cross sections for heavy 
lepton pairs (fi + fi~, or r + r~), can be obtained similarly. The electromagnetic production 
of n + pairs using hadro n be ams was first observed in 63 GeV pp collisions at CERNs 
Intersecting Storage Rings jl07 |. 

For meson pairs like nir, neglecting internal substructure, as is done for Equation 12*61 mav 
be appropriate near threshold. However, at higher pair masses, the quark substructure of 
the mesons be come s important, and the cross section for 77 — > 7r°7r° becomes comparable to 
that for 7c + 7r~ 10^ . In fact, studies of 77 production of mesons pairs are interesting probes 
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of meson structure. Baryon-antibaryon pairs are also of interest, because the reaction probes 
the baryon internal structure. 

Because the cross sections depend on the inverse of the square of the particle mass, 
production of heavier pairs t + t~) is much smaller than for e + e~. Their Compton 

wavelength, Aj = h/rriiC is smaller than the nuclear radius R. This requires the replacement 
L — > £ = ln(7i72<5 / 'rriicR) in Equation l2*6l where 5 = 0.681... is a number related to Eulers 
constant. 

Bound particles, such as positronium or qq mesons are also produced in two-photon 
interactions. The cross section is given by Equation 2. The cross section for 7172 — > X 
depends on the particle's decay width to two photons, T 77 10g| . Since decay and 77 
production use the same matrix el emen ts, only the phase-space factors and polarization 
summations are distinct. One finds 



a{u Xt uj 2 ) = 8n 2 (2J + l)-f^L 8{Au x u 2 - M 2 c 4 ' 



(27) 



where J, M, and TM-ryy are the spin, mass and two-photon decay width of the meson. The 
delta-function imposes energy conservation. 

Using Equation EH the production of mesons with mass M in HI colliders is jsj: 



128 2 ^r 77 

° = — {ZlZ > a) 



[£» + ...] 



(28) 



This equation is obtained by using Equation |21 and the high e nerg y limit (7 ^> 1) of the 
equivalent photon number n{ui) (for more details, see reference HQJ]). 

A more detailed account of the space geometry of the two-photon collision is necessary 
especially for heavier mesons, and will be discussed in Section 4.4. Since spin 1 particles 
cannot couple to two real photons only spin and spin 2 particles should be produced. 

The treatment of bound states in quantum field theory (QFT) is a very complex subject 
(for reviews, see 
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113]). In the case of positronium production by two photons (para- 



positronium) and by three photons (ortho-positronium), standard QFT tec hniques a llow a 
simple and accurate way to calculate the cross sections from first principles |^y|^3]. The 
para-positronium production cross sections a re qu ite large, 19.4 mb and 116 mb, for RHIC 
(Au+Au) and LHC (Pb+Pb), respectively (llj. However, Cou 



these values by as much as 43% for RHIC and 27% for LHC 114j . The cross section for 



omb corrections reduce 



the production of ortho-positronium, which requires three-photon exchange, are also large: 
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11.2 mb and 35 mb, for RHIC and LHC, respectively 114]. Even the ortho-positronium 
cross sections correspond to production rates of 4 and 35 per second respectively. If the 
positronium could be extracted from the the interaction points, they could be used to test 
interesting properties of QFT for bo und states. Relativistic positronium has an unusual 



large transparency in thin layers (see 



llq and references therein). 



The same diagrams for the calculation of the positronium apply for production of bound 
qq pairs ( meson s) in UPC 115j ]. However, proper accounting for the color degrees of freedom 



is needed 
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B. Production of free e + e pairs 

Due to experimental difficulties, Equation |2H1 (and its newer counterparts) has never 
been fully tested. With the construction of RHIC and the LHC, interest in this process has 
grown. For heavy ions, the e + e~ production probabilities are close to one and lowest-order 
perturbative calculations of the cross sections violate unit arity (i.e . d 2 a/d 2 b > 1) 

This observation lead to more detailed calculations E H H, H Q involving 
high-order processes, such as the exchange of multiple photons (Coulomb distortion) and the 
production of multiple pairs, as shown in Figures ^e,f). These processes are important for 
collisions at small impact parameters. Diverse theoretical methods have been considered. 
Perturbative calculations are simple to write down, but they involve rather complicated 
integrals, especially for low energy electrons, due to Coulomb distortion and relativistic 
effects on the continuum electronic wavefunction |9|. A general sum of the contribution of 
diagrams like those in Figure ^e,f) and unitar ity c orrections (involving the production of 
virtual e + e~ pairs) was obtained in Reference [l23j . To account for Coulomb distortions, 
one needs to add to eq. |2H]a term of the form (see eq. 7.3.10 of Reference jsj), 

°c = [f(Zxa) + f(Z 2 a)} a L 2 , (29) 
y7r 



where 



is the Bethe-Maximon correctio n. Eq . EH] was derived in Ref. |l25| . and later confirmed by 



f( x ) = x 2 V — - (30) 

v ' ^ n(n 2 + x 2 ) 

n=l v ' 



independent calculation in Ref. |l2fi| . 



27 



For Pb+Pb collisions at LHC the Coulomb distortion correction reduces the pair- 
production cross section by 14%. Other unitarity corrections further reduce the cross sections 



122 1 - 



by 3% 

The calculation of the production of multiple pairs, as shown in Figure de,f), is directly 
connected with the unitarity problem. It is possible to interpret d 2 a/d 2 b as the mean number 
of pairs produced a t a g iven impact parameter. For Ca-Ca collisions at the LHC (Za = 0.15), 
123], or about 27 000 e + e~e + e~ events per second. In the literature, 



0-2- 



pairs 



0.11 b 



one finds different methods to calculate the cross section for n > 2 pairs. The result of 



Ref. 



1241 ] is a simple fit to numerical calculations. Ref. 1181 ] is based on the probability 



P(b) of the pair production taken from Ref. Ref. |l23j claims that this expression is 



wrong. They derived expressions for this probabi 
123l obtains a n oc L n , whereas References 



ity usi ng tw o different methods. Reference 



1181 ] and |l24j | obtain a n oc L 3n respectively. 



The calculations differ in the method used to include Coulomb and unitarity corrections. 



The production of multip 



approaches. [1J, 112 



125, 



e p airs 



127, 



r as been s tudied w ith a variety of different theoretical 



133, 



134, 



135, 



136 



l3 



The calculation of multiple photon exchange can be considerably simplified in the ultra- 
relativistic limit. In this limit, the electromagnetic field of the ions is squeezed in the plane 
perpendicular to its trajectory (i.e. it can be approximated by a delta-function along this 
plane). In the appropriate gauge, the Coulomb p otential is two -dimensional and the time- 



130, 



138, 



This should be equivalent 



dependent Dirac equation may be solved exactly 
to an all-orders perturbative calculation. 

This approach yields good results as long as ub/'-fhc < 0.1 Above this value the 

delta-function approximation breaks down. Since the most important impact parameters 
for this process are of the order of b ~ h/m e c |9], the calculation can be separated into two 
regions: (a) b ~ h/m e c where the approximation is valid, and (b) large impact parameters, 



for which perturbative calculations are accurate 



130, 



l~T 



138. Il39^ . This method describes well 
the differential cross sections for e + e~ pair production up to energies of the order of 0.l7m e c 2 , 
abov e which the delta-function approximation breaks down for the same reason as above 



140]. The initial calculations using this technique foun d results t hat matched the lowest 
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1321 ] . This was inconsistent 



order perturbation theory without Coulomb correct ions 

with both theoretical expectations and with data jl32j |. However, regularization of the 
integrals was critical; with regularization, the Dirac approach reproduced the lowest order 
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result, with Coulomb corrections 126l Il4l| . This technique allow s for the calculation of 



cross section for free e + e pair production to all orders in Za 141 1 



3 air production has been studied at RHIC in combination with mutual 



Electron-positron 

Coulomb excitation 142 1. As will be discussed in the next section, the mutual Coulomb 
excitations were required to trigger on these events, and also had the effect of selecting events 
with b < 30 fm, where non-perturbative effects were strongest. The cross section, pair mass, 
rapidity and pair v t dis tributions were all in accord with the predictions of lowest order 
perturbation theory jl3(| . The pair p T distribution deviated from the Weiszaecker- Williams 
virtual photon prediction, showing that the photon mass was important in that kinematic 
regime. 

The STAR study suffered from low statistics. Earlier experiments on e + e~ production 
in sulfur-platinum (and sulfur-nuclei) collisions at the CERN SPS had higher statistics, 
)ut l ower beam energies; they also found good agreement with lowest order calculations 
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I 1 

It remains disappointing that these 70-year-old QED calculations are still not fully tested. 
Although many aspects of QED have been tested to high precision, studies involving strong 
fields are much less advanced. Pair-production with relativistic heavy ions (with Za ~ 1) 
is one important example. 



C. Pair production with capture and antihydrogen 



An important phenomenon occurs when the electron is captured in an atomic orbit of the 



n 



projectile, or of th e target |9(. At RHIC and the LHC, this is an important source of beam 



particle loss llOj]. The produced beam of single-electron ions carries considerable power 
14f| : at the LHC, at full luminosity, the produced +81 Pb beam carries sufficient power 
to quench the s uperconducting accelerator magnets; this limits the LHC luminosity with 



heavy-ions 



146]. 



One striking application of this process was the recent production of antihydrogen atoms 
using relativistic antiproton beams ■?]. The positron is produced and captured in an orbit 
of the antiproton. The expression 

ZfZla 8 H 2 1 



a 



3.37T- 



exp (27rZ 2 a) 



(L — 2.051) 



(31) 
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for pair production with electron capture in the nucleus with charge Z2 is obtained in first 
order perturbation theory j^. Although Equation |^ works reason ably well for explaining 
antihydrogen production, it is only valid for small Z^ (Z t < 15) [7], 147 . Il4s| |. For large Zj, 



as w i th the experiments at R HIC and LHC, non-perturbative calculations may be necessary 

152| . EquationEUincludes higher order effects related to the electron 



119 



130 



149, 150 



151 



capture, but is not a complete all-orders result. The additional higher order corrections are 
apparently small, and Equation should be usable for most purposes. The fraction with 
the exponential term is due to the distortion of the positron wavefunction. It accounts for 
the reduction of the magnitude of the positron (continuum) wavefunction near the nucleus 
where the electron is localized (bound). 

Equation I3*T1 shows that the cross section depends on energy as a = A In 7172 + B, where 
the coefficients A and B depend on the nuclear charges. This scaling was confirmed in 
numerical cal culat ions of Reference jl3o| and was used in the analysis of the experiment 
in Reference Jl53l . shown in Figure IHl^a). The comparison between theory and the data 
of Reference |l53| . is not completely valid since atomic screening was not included. When 
screening is present the cross sections are smaller by at least a factor 2-4 (see Equation 7.4.3 
of Reference 



0; 



Similar reactions occur with hadron pairs. Processes like pion pair production with cap- 
ture (or, similarly, photoexcitation of a A resonance, which decays by 7r + emission) increase 
the nuclear charge by 1 and can turn lead into bismuth, plus a tt + . This change cannot 
occur electromagnetically, since e + will not bind to lead. The se transmutations have been 



studied at the SPS/CERN for 208 Pb at 160 GeV/nucleon 
quantitatively with electromagnetic excitation calculations 



154 

i [87 



. Th e data can be described 
Il55| . For high-Z nuclei, the 



dominant contribution to nuclear-charge pickup is due to electromagnetic production of tt~ 
by virtual photons. The electromagnetic contribu tion is completely negligible for a similar 
experiment at an energy of 10.6 GeV per nucleon jl5fij |. 



D. Two-photon production of mesons 

For the production of composite particles, one can use Equation |2H1 as a first guess. 
However, even the lightest mesons (tt°) require photons of relatively large energy (> 70 
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FIG. 8: (a) Pair production with capture for Pb 82 + (33 TeV) beams on several targets. The solid 



and dashed curves are theoretical calculations 



1531 ] . See text for more details, (b) The tagged total 



photon flux (accompanied by a single giant dipole resonance excitation) for a complete calculation 
(solid line) and a simplified 'box' integration (dashed line), compared to the untagged flux (dotted 
line), for gold at RHIC. Although tagging reduces the low-energ y fl ux by an order of magnitude, 
at high energies, the difference is much smaller. From Reference |22(. 

MeV). Mesons are produced primarily in collisions with relatively small impact parameters 
(compared to 2Ra) where hard photons are more abundan t. Su bstantial changes are required 
to Equation |2H to account for the collision geometry jrsl . Il5?| . 
One can rewrite Equation |2] more conveniently as 



ox = I ds^^- *^x(W yy ) , (32) 



where = Aujiuj 2 is the square of the center-of-mass energy of the two photons, 

a^^xiWyy) is the two-photon production of particle X, and dC/dW J7 is the "photon- 
photon luminosity". dC/dW^ can be multiplied by the ion beam luminosities, yielding an 
"effective" two-photon luminosity dC e ff/dW J7 which can be directly compared to other two- 
photon luminosities, such as at e + e~ or pp colliders 158| . With heavy ion beams, the LHC 
two-photon luminosities are much higher than are available elsewhere, either with proton 



beams at the LHC or at the e + e LEP-II collider for energies up to \fW^ ~ 500 GeV (14 1. 
Table ITVl shows the cross sections for the production of C = ev en mesons f or the 



RHIC and LHC colliders. Other calculations were done in References |8£ 



100, 153. The 
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TABLE IV: Cross sect ions for two-photon production of (C = even) mesons at RHIC (Au+Au) 



and at LHC (Pb+Pb) 



115]. 



meson 


mass [MeV] 


a RHIC j mb j 


a LHC [ mb ] 


7Tn 


134 


4.9 


28 


Ti 
i 


547 


1.0 


16 


rf 


958 


0.75 


21 


/ 2 (1270) 


1275 


0.54 


22 


a 2 (1320) 


1318 


0.19 


8.2 


Vc 


2981 


3.3xl0~ 3 


0.61 


XOc 


3415 


0.63xl0~ 3 


0.16 


X2c 


3556 


0.59xl0~ 3 


0.15 



cross section corresponds to a 77 — > tt° rate of 30 events/second with lead beams at the 
LHC. For heavier mesons, like i] c , the rate is still large, of the order of 1 per minute. 

For mesons of comparable mass, the two-photon cross sections in Table HVl are about two 
orders of magnitude lower than the cross sections for photonuclear vector meson production 
(Table |H]) . This difference stems from the different coupling strengths of the strong and 
electromagnetic interactions, a s ~ 1 and a em ~ 1/137. 

Two-photon meson spectroscopy is thus greatly complicated by the large background 
from photonuclear interactions. For example, with lead beams at the LHC, the rate of J/ip 
photoproduction followed by J/tp — > ^r\ c is about 2.5 per minute, higher than the 77 — > 7] c 
rate. 

Although it may be possible to separate the different event classes with cuts on meson 
Pt, rapidity, and final state particles, the vector meson background seems daunting to most 



efforts 
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E. Searches for New Physics 



The LHC will reach high enough energies that two-photon interactions will be an attrac- 
tive place to search for some types of new physics. Many early calculations focused on the 



search for the Higgs |163l . 11641 Il65j . Other examples include super symmetric particle pairs, 
magnetic monopoles and possible extra spatial dimensions 1J]. The LHC will also be able 
to probe vector boson couplings through reactions like 77 — ► W + W~ . 

The two-photon production rate for the Higgs is small enough that, for most models, it is 
likely to be discovered in hadronic interactions. However, for standard model Higgs masses 
under « 200 GeV, with medium ion beams, the 77 channel should produce a handful of 
events per year [14|. In some sup ersym metric scenarios, the production of the Higgs in UPCs 
could be significantly enhanced |l59|. The 77 production channel could also be studied in 
pp collisions; the greatly increased luminosity and running time will more than compensate 
for the smaller production cross section. However, in pp collisions, there is a considerable 
background due to diffractive interactions. It may be possible to separate 77 from diffractive 
interactions by studying the px of the scattered protons. This could be done by placing small 
detectors, kn own as Roman pots inside the beampipe, to detect protons scattered at very 



small angles 



160]. 



Supersymmetric particle pairs are similar story. If present, they are likely to be discov- 
ered in hadronic interactions. However, if supersymmetry is correct, a large number of new 
particles are likely to be present, and the 77 production of sparticle pairs is likely to pro- 
vide significant new information; two photon production is sensitive to significant regions 



of phase space 



159]. Two-photon interactions are sensitive to sparticles that do not partic- 



ipate in the strong interaction. Photonuclear interactions may also be useful for studying 
supersymmetry. 

Real or virtual magnetic monopoles can be produced by two-photon fusion. The CDF 
collaboration searched for the process 77 — ► 77 wat the Fermilab Tevatron, and set mass 
limits on magnetic monopoles with various charges |l6l| . The LHC will be able to do far 
better. 

The presence of extra dimensions could be detected via the two-photon production of 
gravitons. The cross s ectio n to produce a graviton, 77 — > G increases in the presence 
of compact dimensions |l62| . There are unresolved theoretical issues regarding the cross 
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section, but one calculation finds that, for 2 extra compact di mens ions, at the LHC, the 



cross section is of order 1 nb for lead and 10 pb for calcium 162]. Rates are not given 
for proton beams, but, for calcium, and likely protons, a few events would be produced 
each month. The experimental signature of graviton production has not been worked out in 
detail. 

Most of these 'new physics' channels involve relatively high px particles, and so should 
be within the purview of the planned trigger setup for the ATLAS and CMS detectors. This 
may not be true for supersymmetric final states; two charged sleptons that don't interact 
hadronically will challenge any trigger. 

V. MULTIPLE INTERACTIONS BETWEEN A SINGLE ION PAIR 

Because heavy-ions have such large charges, a single ion pair can undergo multiple electro- 
magnetic reactions in a single interaction. Even though the reactions may be independent, 
the geometry introduces correlations between the photon energies and polarizations. Mul- 
tiple interactions are also a key experimental tool, allowing for many cross checks under 
different triggering conditions. One example of such a reaction is the photoproduction of a 
p° meson, accompanied by the mutual Coulomb excitation of both nuclei: 

Au + Au -> Au* + Au* + p°. (33) 



This reaction was studied by the STAR collaboration [251 ] . This process occurs predomi- 
nantly via 3-photon exchange, as is shown in Figure Q^g). 

The STAR collaboration has also observed four photon react ions , such as the production 



of an e + e pair accompanied by mutual Coulomb excitation 1421 ] . as is shown in Figure 

raw- 

In a multi-photon process, each photon emission may be treated independently, if the 
energy lost by the nucleus is not significant. As long as the photon emission does not excite 
the emitter, the reactions may be treated as completely independent. The cross section is 
calculated in impact parameter space 

a = I d 2 bP(b) (34) 
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where P(b) is the probability for the reaction to occur at impact parameter b. This is 

P(b) = [ —N{u,b)a lA {uj). (35) 
J u 

When the cross section for a reaction is very large (as with e + e~ production or GDR ex- 
citation) the naive P(b) calculated in Equation 1351 mav exceed 1. P(b) should then not be 
interpreted as a probability but rather as the mean number of produced particles at that 
impact parameter. The actual number of produced particles follows a Poisson distribution 
with this mean. The generalization to 2 (or more) photon exchanges is obvious. Calculations 
using this approach accurately predicted the cross section and kinematic distributions for 
Reaction 1221 25]. 

This factorization only holds if several conditions are satisfied. Photon emission must 
not excite the emitting nucleus and the photons must be emitt ed in dependently. As long as 
the fractional energy loss of the nuclei are small, this is valid |l68^ . Finally, the excitation 
must not change nuclear form factor significantly on the time scale of the reaction (the 
'frozen nucleus' approximation). As long as these strictures are satisfied, the ordering of 
the subprocesses is unimportant. These conditions hold for heavy-ion collisions. For proton 
beams, with rj = Z\Z2e 2 a/ (3 <C 1, the factorization is on weaker ground, because of the 
poorly defined impact parameter [22]. 

It can be convenient to treat one reaction as a trigger (or selector) for a range of impact 
parameters. Picking events with mutual Coulomb excitation, for example, selects events 
with small impact parameters. The reason can be seen in Equation |^ In the low-energy 
limit (u <C jHc/Ra), for a fixed u, P(b) ~ N(u,b) ~ 1/b 2 . For a two-photon reaction, 
Pl{b)p2{b) ~ \/b A . The mean impact parameter b n for an n-photon interaction is [3] 

/d 2 6ftP 1 (6)...P w (6) 
bn ~ J<PbP 1 (b)...P n (b)- W 

Here, b m i n = 2Ra is the minimum impact parameter, and b max = ^K/Ra- For n = 1 the 
result b = (b max — b min )/ In (b max /b min ) is not so useful. However, for 2 or more photon 
exchanges 

b n >\ = ^Ra- -', (37) 

2n — 3 

bmax drops out, leaving a simple result. For n = 2, this reduces to 6 2 = 2Ra; for larger 
n, b is even smaller. At heavy-ion colliders, mutual Coulomb excitation is an effective 
trigger for selecting low-impact parameter events. Detailed calculations of the median impact 
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parameter in 1 and 3 photon interactions find a similar scaling [ 76] . Reducing b is very helpful 
in studying interference in vector meson production, by increasing the p? range over which 
the interfence is visible. 

This selection can also be viewed in momentum-space. In the low-energy limit, the 
photon flux n(u) (Equation^ scales as 1/u. However, when an additional photon is present, 
the spectrum becomes much harder. The spectrum for photons that are accompanied by 
Coulomb excitation is: 

pfic/u Const 
n(u)m \ d 2 bN(uj,b) — -— w Const. ; (38) 

J2R A b 2 

The extra photon adds a 1/b 2 weighting, and the resultant flux is independent of the photon 
energy. 

Figure Efb) compares the spectra with and without tagging. By selecting reactions with 
additional accompanying photon interactions, experimenters can 'tune' their photon beam, 
hardening or softening the spectrum. This 'tuning' allows many cross-checks. For example, 
in vector meson production, y = In (2u;/Myc 2 ), but there is a two- fold ambiguity over which 
nucleus emitted the photon. By comparing vector meson production with and without mu- 
tual Coulomb excitation, it is possible to account for this ambiguity and find the production 
cross section function of photon energy. 

The coupling is also very useful in experimental triggering. A simple reaction like mul- 
tiple Coulomb excitation can be used to trigger on small-impact-parameter collisions; these 
remainder of the event can then be studied without experimental bias. 

For two-photon final states, like e + e~ pairs, the situation is more complex, because the 
particles are produced outside the nuclei. However, two-photon reactions are also enhanced 
at small nucleus-nucleus impact parameters. The STAR collaboration used mutual Coulomb 
excitation to study low-mass e + e~ pair production at RHIC, since it was not possible to 



trigger on the e + e~ pair itself [1421 ] . The presence of the mutual excitation also significantly 
hardens the pair mass spectrum. 

In multiple photonuclear interactions, the photon polarizations are collinear. Photons 
are linearly polarized along the electric field of the emitting nucleus. In photonuclear inter- 
actions, the electric field vector follows the impact parameter vector. When a single nucleus 
emits multiple photons, these photons all have the same linear polarization. When the other 
nucleus emits a photon, it will have the opposite polarization. 
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When the final states are sensitive to the photon polarization, then angular correlations 
will be present. In p° decay, in the plane perpendicular to the p° direction, the angle 
between the photon polarization vector and the n + (or 7r~) direction is distributed following 
a cos 9 distribution. For the case of two independent p° decays, with uncorrelated photon 
polarization, the angle between the two n + , 5<p is evenly distributed between and 2n. With 
the polarization correlation, the angular correlation is 

N(A(f)) cos(2A0). (39) 

A similar distribution is also expected for neutrons from giant dipole resonance. Correla- 
tions between neutron in p° production and GDR excitation(s) should also be measurable 
at RHIC. In the longer term, GDR excitation neutrons could be used to tag photons ac- 
cording to their polarization direction, allowing for studies with polarized photons. Similar 
polarization should also occur for medium-energy nuclear reactions. 

In addition to the correlations due to the geometry, multiple interactions may be a place 
to study Bose-Einstein correlations, such as in p°p° production by two independent pho- 
tonuclear interactions. When the two p° are produced on the same nucleus, the production 
should be bosonically enhanced (exhibit super-radiance) when the p have (in the nuclear 
target rest frame) a momentum difference Ap < H/Ra- 

VI. EXPERIMENTAL POSSIBILITIES AND LIMITATIONS FOR ULTRA- 
PERIPHERAL COLLISIONS 

The experimental study of electromagnetic interactions at high energy colliders is quite 
new. Since the characteristics of these interactions are very different from the more common 
hadronic interactions, most existing and planned detector systems are not optimized to 
study them. This section will discuss some of the general experimental possibilities and 
limitations at current and future colliders. 

At hadron colliders (in contrast to electron accelerators) it is in general not possible to 
detect the outgoing projectiles following an electromagnetic interaction. This is because of 
the small momentum transfers involved. Tagging of nuclei will never be possible because the 
angular deflection following the coherent emission of a photon is smaller than the angular 
dispersion of the beam. Proton tagging has been proposed, but requires extremely high 
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FIG. 9: Transverse momentum distribution for events with exactly 2 reconstructed charged pions. 
Data from the STAR collaboration [7^. The Monte Carlo calculation is based on [ 7^ . 



resolution Roman pots 



160] 



For exclusive particle production, characterized by the emission of only a few final state 
particles, good signal to background ratios can be achieved by selecting events with small 
total transverse momentum when the event is reconstructed [88(. The total event transverse 
momentum is the sum of the momentum transfer from each projectile, which is determined 
by the form factors and can be calculated accurately. The method works best for heavy 
nuclei. This is illustrated in Figure EJ which shows the px distribution for exclusive p° 
production in Au+Au interactions at RHIC. The coherent peak for events with exactly one 
reconstructed tt + and one reconstructed 7r~ can be clearly seen with < 100 MeV/c. The 
incoherent background can be estimated by measuring events with two reconstructed pions 
with the same charge (dashed histogram in the figure). 



The technique of using the total event pt to identify electromagnetic interactions works 
only for coherent and exclusive particle production. Other photonuclear interactions require 
different approaches. Incoherent photon-induced interactions are characterized by a gap in 
rapidity void of particles between the rapidity of the photon emitting projectile and the 
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rapidity of the produced state. This distinguishes the electromagnetic interactions from 
ordinary hadronic interactions. 

For hadronic interactions, the probability of having a rapidity gap of width Ay void of 
particles decreases exponentially with increasing size of the interval: 



where (dn/dy) is the mean multiplicity per unit of rapidity. Requiring a rapidity gap of 
width Ay = 2, leads to a hadronic rejection factor of 10~ 2 to 1CT 3 for single nucleon-nucleon 
interactions at RHIC and LHC energies. Additional rejection of hadronic events can be 
obtained if the fragmentation of the beam nuclei is detected, for example in forward ( "Zero- 
degree") calorimeters. In hadronic interactions both nuclei normally break up, whereas in 
photoproduction the photon emitter normally remains intact. Further discussions of the 
experimental aspects can be found in [^j]. The conclusion is that sufficient rejection can 
be obtained with current and planned detectors at heavy ion colliders to study, e.g., the 
production of heavy quark pairs. 

One significant background to electromagnetic diffraction is hadronic diffraction. 
Pomeron exchange can also produce rapidity gaps, and double-Pomeron inte racti ons pro- 
duce isolated systems in a central detector, with rapidity gaps on each side jlOo| . These 
backgrounds are not an issue for vector meson production, because double-Pomeron inter- 
actions do not produce J PC = 1 final states, but these backgrounds could be problematic 
for other reactions. Because the strong force has a short range (< 1 fm), Pomeron exchange 
between nuclei and Double-Pomeron interactions can only occur between surface nucleons 
in grazing collisions; this limits the cross sections. Unfortunately, quantitative estimates 
are lacking. However, for pp collisions, Pomeron interactions may constitute a significant 
background. 

Another experimental challenge is finding an efficient trigger with adequate background 
rejection. Since the outgoing beams are not tagged, it is necessary to trigger on the particles 
emitted from the final state X to study the reaction A+A — > A+A+X, . These particles 
are usually produced at or near mid-rapidity. Most of the current and planned heavy- 
ion experiments at RHIC and the LHC have triggers that are optimized for hard, high px 
interactions. These triggers are difficult to adapt to low multiplicity low-pr final states. 
STAR at RHIC is an exception, but it faces challenges because of the low allowable trigger 




(40) 
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rates [169(. Backgrounds from beam-gas interactions, grazing hadronic collisions, ambient 
neutrons and other beam-related backgrounds are serious concerns to any experiment. 

In combination with the large probability for multiple Coulomb interactions with heavy 
beams, factorization makes triggering on neutrons emitted in the forward direction following 
single or double Coulomb excitation an attractive alternative. This reduces the photon flux, 
but provides some control of the impact parameter distribution, as discussed in Section 5. 

The experimental feasibility of studying ultra-peripheral collisions is demonstrated 
through the data presented earlier. However, the considerations above limit the types of 
processes that can be studied. The study of two-photon production of multiple e + e~-pairs 
seems difficult because of the extremely low px of the emitted electrons. A dedicated ex- 
periment would do better, but that is likely to be a difficult proposition at a high-energy 
physics laboratory. The study of meson spectroscopy at the LHC also seems problematic 
because of the lack of triggers for low momentum charged particles around mid-rapidity. 



VII. CONCLUSIONS 



We have presented the formalisms for studying photoproduction and two-photon reac- 
tions at hadron colliders, and discussed some of the more interesting applications of these 
techniques. Low-energy nuclear physics has used these techniques for many years, with good 
results. However, only recently have higher-energy machines like RHIC and the Tevatron 
begun to study particle production in very strong fields. 

The amount of experimental data on UPCs is still rather limited. Despite this, the field 
has developed tremendously over the past decade, with much theoretical progress. The new 
data from RHIC is helping to focus new theoretical work on the channels that are most 
readily accessible. With the new results from RHIC and the coming LHC startup, UPCs 
are now ready to make substantial contributions to many areas of physics. 

The LHC will produce photonuclear interactions with 10 times the energy available at 
other accelerators like HERA. This will open a huge window to search for new physics 
processes (some of them not accessible in hadronic collisions), measure gluon densities at 
very low Feynman-x, and perform a host of other measurements. At the same time, UPC 
reactions are very important for machine operations - e + e~ production with e~ capture will 
limit the LHC luminosity with heavy ions. 
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